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Abstract 

A crucial issue in relativistic plasma, particularly relevant in the astrophysical context, is the 
description of highly magnetized plasmas based on a covariant formulation of gyrokinetic dynamics. 
An interesting case in question is that in which the background electric field (produced either by 
the same plasma of by other sources) results suitably small (or vanishing) with respect to the 
magnetic field, while at the same time short-wavelength EM perturbations can be present. The 
purpose of this work is to extend the relativistic gyrokinetic theory developed by Beklemishev et 
al. [1999-2005] to include, in particular, also the treatment of such a case. We intend to show 
that this requires the development of a perturbative expansion involving simultaneously both the 
particle 4-position vector and the corresponding 4-velocity vector. For this purpose a synchronous 
form of the relativistic Hamilton variational principle is adopted. 

PACS numbers: 52.30.Bt, 52.55.Hc, 52.55.Dy 
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I. 1. INTRODUCTION 



This work is a part of a systematic investigation concerning the relativistic gyrokinetic 
theory developed by Beklemishev et al. and its apphcations. In this paper, in 

particular, we intend to propose an extended formulation of covariant gyrokinetic theory 
specifically intended to include the effect of 4-velocity perturbations, as appropriate for 
the description of relativistic magnetoplasmas in the presence of suitably "small" ( and 
regular, but otherwise arbitrary) electromagnetic and gravitational perturbations produced, 
via collective interactions, by the same plasma. The physical motivations of the investigation 
are related to the occurrence of strongly magnetized relativistic plasmas in curved space- 
time. In the astrophysical context these plasmas are typically related to the existence of 
accretion disks, plasma inflows and outflows and relativistic jets, all occurring typically 
close to massive object, such as neutron stars, black holes (BH) and active galactic nuclei 
(Mohanty ^). These plasmas can be characterized by the presence of ultra-strong magnetic 
fields, possibly generated by the same plasma. Although, direct measurements of B can be 
very difficult, for example in jets there is yet no reliable observation of the magnetic field 
intensity, the partial information already available for other cases indicates that in suitable 
circumstances the magnetic field can become extremely intense. For example, sufficiently 
close to the surface of pulsars the magnetic field intensity can reach values well above the 
so-called quantum critical field strength of B^. = mlc^/eh = 4.4xlO^^G. Magnetic fields of 
this order or larger are said to be super strong. A similar situation might occur for BH's 
surrounded by strong magnetic fields, and in particular in the regions of space-time suitably 
close to the so-called photon-sphere [lO|, where coUimated outflows can arise. 

While such fields are typically expected to manifest together with intense electric fields 
(E), there can be cases in which there results in particular 

E/B-^l (1) 

{weak electric field assumption) . Here E and B are the invariant eigenvalues of the elec- 
tromagnetic (EM) field tensor F^^ (see below). This may happen, for example, in plasmas 
in which the magnetic field is produced primarily by strong diamagnetic currents. These, 
by definition, are electric currents flowing in the plasma which are driven by gradients of 
the relevant fluid fields (in particular, fluid velocity, number density, temperature, magnetic 
field intensity, etc.) and, instead, are independent of the (local) electric field. 



An example is provided by plasma outflows (or jets) in which the kinetic energy of the 
escaping particles is usually believed to by provided by the presence of a quasi-stationary 
magnetic field moving with the plasma [isl- In all such phenomena EM turbulence is 
expected to be ubiquitous and to affect in a significant way plasma dynamics. The region 
where these processes may act is typically too close to the central object (the BH or a 
neutron star) to be resolved observationally. Thus, one possible way to obtain some insight 
into this region is by studying the properties of these phenomena. In particular, the emitted 
radiation spectrum could be significantly modified by the effect of a finite curvature of space- 
time, as well as by short-wavelength EM and gravitational perturbations in the region of 
emission fl^. 

An important theoretical problem is therefore the prediction of such effects based on 
single-particle gyrokinetic dynamics. Such phenomena require their description in the frame- 
work of a consistent covariant theory taking into account the effect of space-time non- 
uniformity as well as the possible presence of wave-fields, in principle, both EM and gravi- 
tational. A convenient formulation is provided by the covariant gyrokinetic theory recently 
developed by Beklemishev et al. [ll, Q]- As well known, object of particle gyrokinetics is 
the description of particle dynamics in terms of a reduced set of variables (so-called gyroki- 
netic). These include, in particular, the fast gyration angle around the magnetic flux lines, 
i.e., the gyrophase 0, which results ignorable. Gyrokinetic theory is therefore realized by 
the construction of a suitable, generally non-canonical, phase-space transformation. From a 
mathematical viewpoint, gyrokinetic theory provides a perturbative description of particle 
dynamics in terms of "simpler" equations of motion, which are characterized by a reduced 
number of effective degrees of freedom. Goal of the present paper is to extend the Bek- 
lemishev and Tessarotto formulation to obtain the gyrokinetic theory in curved background 
space-time which applies also to the special case [1] and includes the simultaneous effect of 
short-wavelength EM and gravitational perturbations. Starting point is the Larmor-radius 
ordering, requiring e ~ -p^ ~ -C l,and A = r^!<^ = ^ ^ 1, where rr is the Larmor 
radius and L e m, Lg, lespectivelj denote the characteristic scale lengths of the background 
and wave EM and gravitational fields. In particular, we intend to determine new expressions 
for the magnetic moment = Ji and of the cyclotron frequency ^0, possibly modified by 
the effect of space-time curvature. 
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II. 2. THE SYNCHRONOUS HYBRID HAMILTON VARIATIONAL PRINCIPLE 
IN CURVED SPACE 



The starting point for the construction of the gyrokinetic theory in curved space-time 
is the proper treatment of the variational principle which determines the gyrokinetic world 
lines {r'"(s),s G /} in the presence of a non-uniform metric tensor (7^,y(r°),as well as the 
EM 4- vector potential A^. Here r"^ denotes the gyrocenter position 4- vector (in the sequel 
the apex will denote gyrokinetic variables or phase-functions evaluated at the gyrocenter 
position). Historically (see for example Einstein [sl, Landau & Lifschitz 6]) the relativis- 
tic equations of motion were obtained using an asynchronous variational principle, stem- 
ming from the minimum-action principle (relativistic Hamilton principle). Accordingly, 
the trajectory of a point-particle is described as an extremal curve (i.e., a geodesies) in 
the 4-dimensional space-time. The choice of an asynchronous formulation, rather than an 
asynchronous one, is a natural consequence of the fact that the variational principle is a 
generalization of the well-known Hertz-Fermat-Maupertius variational principle (Goldstein 
et al. 81]) and that the 4-velocity of a particle, u^{s) = ■^r'^is) , is a unit 4-vector, i.e., it 
must fulfill the constraint u^u'^ = 1 {realizability condition) , where = g^uu'^- As a conse- 
quence, it follows that the line element ds along the world line {r^{s),s G /} must satisfy 
the constraint ds"^ = gf^^dr'^{s)dr''{s), which implies that the asynchronous variation of ds 
results A {ds) 7^ 0, A denoting the operator of asynchronous variation. On the other hand, 
adopting the [relativistic] asynchronous principle one has to face the conceptual difficulty 
represented by the fact that - contrary to what happens for the non-relativistic Hamilton 
principle - it cannot provide a variational principle for the relativistic Hamilton equations 
(at least using a Legendre transformation like in non-relativistic mechanics). In reality, the 
difficulty can be circumvented (see for example Wheeler et a/.j3] and Goldstein et al. ^) 
adopting a suitable synchronous variational principle. In such a case it is immediate to show 
that the canonic form of the relativistic Hamilton equations is recovered, just like in non- 
relativistic mechanics, by means of a suitable Legendre transformation. Finally, it must be 
recalled that variational principles, both synchronous and asynchronous, can be expressed 
in arbitrary "hybrid" (i.e., generally non-canonical) variables, provided these variables are 
connected by a diffeomorphism. to the original one [i.e., for example to be identified with 
the canonical variables], thus yielding a hybrid variational principle. These variable can 
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also be superabundant, since the hybrid variables can be defined in such a way to span, in 
principle, phase-spaces of arbitrary dimension. Thus, they can be described in the frame- 
work of the so-called constrained dynamics (Dirac 9|), in turn relying on the well-known 
method of Lagrange multipliers. In particular, as a special choice, the hybrid variables can 
also be identified with appropriate gyrokinetic variables. In this section we intend to obtain 
a synchronous variational principle of this type, in which the hybrid variables are provided 
by a suitable set of superabundant non-canonical variables. For definiteness let us introduce 
the action S{y) = 7(y) where y denotes the set of variables y = {r^, u^, A) , 7 a suit- 
able differential 1-form and F(y, s) an arbitrary gauge function [i.e., a real function which 
is at least C*^^-* in its existence domain]. In particular, let us impose that the fundamental 
differential 1-form 7 takes the value 

7(y) = 9f.u [qA" + u"] dr'^ + A [g^^u'^u" - 1] (2) 



Here the notation is standard 



2|. Thus, A'^ are the counter- variant components of the EM 
4-potential and q = Ze/ nioC is the normalized charge, and, moreover, r^ and coincide 
with the countervariant components of the 4- vectors r (4-position) and u (4- velocity). It is 
immediate to show that phase-space trajectory of the point particle with 4-position r" and 
4- velocity results, by construction, an extremal curve of the variational functional S'(y). 
Another equivalent possibility is to represent the differential 1-form in terms of 

7(r^ u,, A) = [qA^ + u^] dr'^ - ^ [g^-'u^u, - 1] , (3) 

where the Lagrange multiplier A has been set equal to its extremal value, A = —1/2. Finally, 
let us require, in case ([2]), that 5* be defined in the synchronous functional class of real 
functions 

{y} ^ {x(s) = {r'^{s),u'^{s),X{x)) : x(s) G C^'\l), (4) 
x(s) eR^,I= [S1,S2] CM,x(s,;) = xi(z = 1,2)}, 

while for ([3]), A is considered prescribed (A = —1/2) It is immediate to prove that in both 
cases and the variational principle 

6S = (5) 

delivers the correct relativistic equations. We remark that the synchronous variations 
6r^, 5 a must be all considered independent. Indeed, the physical realizability condi- 



tion is satisfied only by the extremal curve. Hence, in a proper sense, this represents an 
example of unconstrained variational principle. 



III. 3. EXTENDED GYRODYNAMICS 

To obtain the new extended formulation of gyrokinetic theory, the variational functional 
S'(y) must be expressed in terms of gyrokinetic variables, which requires a suitable extended 
phase space transformation to a new set of variables y\ including in particular the gyrophase 
(f) which describes the fast gyration motion of particles around magnetic flux lines. Usually 
such a transformation is obtained by means of perturbation theory in e, to be considered 
infinitesimal and defined as the ratio between the Larmor radius and the smallest character- 
istic scale length of the magnetoplasma. In particular this requires that the EM 4— potential 
A'^ and the metric tensor ^f^^ are assumed to be expressed in terms of the representations 

= iA^(r-)+ea^(r7e), (6) 
9,. = G^.K,e)+e^7;.(r70 (7) 

{Larmor-radius ordering). Here ^A^{r°',€) and G^uir°',€), denoting the equilibrium terms, 
are assumed to be C^°°^ in all variables, while ^a^{r'^/^) and ig]^y{r'^ /C)i define the short- 
wavelength (and high-frequency) perturbations. In particular, is an infinitesimal parame- 
ter, assumed generally independent of e, and in particular such that e -C A, which charac- 
terizes the tensors ia^{r" H) and ,^(7^^(r^/^). Hence, the equilibrium quantities can be Taylor 
expanded near r'" 

^ = ^ + ^^ + o(e), (8) 
G,, = G'^, + eGl, + o{e), (9) 

where A'^" = (r'°) , G'^^ = g^^ (r'") , A'^ = r'^d'^A''' and Gj^^ = r^d'^G'^^. The perturba- 
tions of relevant fields depend generally on the gyrophase angle through er" (y*) . The per- 
turbative scheme converges, at least in an asymptotic sense, only if one assumes the validity 
of a suitable ordering conditions (Larmor radius ordering), implying that EM contributions 
must dominate, in a suitable sense, over inertial ones in the variational functional. Formally 
this corresponds to replace the charge e by e/e, leaving the rest mass mo unchanged. In 
this work, extending the approach developed by Beklemishev et al. , the gyrokinetic 



transformation is taken of the general form (r",M^) ^ = (r'°,M'°) , where r'" and u'" 
identify respectively the so-called gyrocenter 4-vector and a transformed 4-velocity, to be 
determined respectively in terms of the two power series 

oo 

= r'^' + J^^'r^ (10) 

s=l 

CXD 

= u"^©^r<(y) (11) 

s=l 

[extended gyrokinetic transformation) . Here rf (y*) and (y*) (z = 0, 1,..) are respectively 
the 4-position and 4-velocity perturbations, while the operator © must be suitably defined 
as corresponds to the 4-velocity addition law In particular, the perturbations r^{y'^) 
and vf (y*) are defined in such a way that both are pure oscillatory with respect to the 
gyrophase 0'. As usual, m'" can be characterized in terms of suitable independent variables, 
one of which includes the gyrophase 0'. This is obtained by projecting u"^ along the four 
independent directions defining the so-called fundamental tetrad unit 4-vectors (with 
a = 0, ...,3), hereon also denoted {l''^ , l"'^ , l'^ , t^) . In particular, the latter can always be 
identified with the EM fundamental tetrad, i.e., the basis vectors of the EM field tensor 
admitting , in this case to be evaluated at the gyrocenter position, F'^^^ = F^yi^x'""), with 
eigenvalues B and E. In particular, the space-like vectors 62,6^^3 are assumed to satisfy the 
eigenvalue equations 

= Be,,, (12) 
Fl^el = -Be2,. 

Here = g'^^e'^ are the covariant components of and the labels a assume a tensor 
meaning in the tangent space. Thus we shall raise and lower them in terms of = //"^e;,^, 
by means of the Minkowskian tensor rj = diag (1,-1,— 1,-1). It follows to leading order in 

£ 

u'^" = a^'cos(j)' + b>'sm(f)' + (13) 

which can be taken as a definition for the gyrophase 0. Here 'g = {g)^ = J^^ g (y) dcj)' 
denotes the 0'— average of a function g (y) and a'^ and can be identified, up to a constant 
factor, with the two space-like 4-vectors, 63 and 63 (which together with are evaluated 
at the gyrocenter position x'"), thus letting = w ei^ and — wcg. The components a^, b^ 
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and m'^ are not completely arbitrary, but satisfy the constraint g^^u^u^ = 1 for all 0'. Thus, 
we shall demand the mutual orthogonality 

af'b^ = a^^u^ = Vu^ = 0, (14) 

while their moduli are expressed via a single real scalar w 

a^a^ = = 1 - (7^.M^r = -w^. (15) 

We stress that only to the leading order = g'^Ja^ where g'^^^ = g^^{x'°') (and similarly 
ttfj, = g'^j,a" , = g'^J)^), while in general = g'^^u^ ^ g^vu'^ ■ From orthogonality relations 
f fn|) one finds the countervariant representation: 

M^ = MOe(; + Mlle^. (16) 

Here = uq, and m" is the "parallel" component m" = —u\\, with uq and u\\ denoting 
the covariant components of u'^, and eg, ei are the remaining time-like and space-like unit 
vectors of the tetrad. From Eq. (ITSl) one obtains Mq = 1 + m| + w^, which permits to express 
Mo by means the independent functions u\\,w. With the above positions for the transformed 
4- velocity it can be shown that (y*) are purely oscillatory arbitrary 4- vectors, namely the 
0'— averages of rf (y*) are zero. 



IV. 4. GYROKINETIC EQUATIONS IN CURVED SPACE 

The extended gyrokinetic transformation (fTOl) . (ITTl) can be carried out in principle at any 
order in e and ^. Here, to illustrate some interesting features of the approach we intend 
to formulate the leading-order perturbative theory with respect to the parameter e (i.e., 
of order e^). In this case the variational functional ^(y) when expressed in terms of the 
gyrokinetic variables (^(f)', Uq, u\\,fj,, r'^, can be shown to take the form 

Sg{y') = j 7,(0',Mo,M||,^,r'^A), (17) 



where 



73 (^0',Mo,M||,/l,r''^,A) = |0A;, + n||/^ + Mov) rfr'^+ (18) 



is the fundamental gyrokinetic differential 1-form. In particular, in case ([3]) there results 
A = —1/2 in Eg. ffTSl) . Here the notation is standard and in particular u is the relativistic 

Ti p n 

magnetic moment pj, |2, 13|. Next, we wish to investigate the consequences of the new 
synchronous variational principle introduced in Sec. 2, in which the . We have proven that 
the Lagrange multiplier A must necessarily take the value A = — | in both cases ([2]) and 
([3]), thanks to the variational principle ([5]). The same result manifestly holds also for the 
gyrokinetic functional ^^(y'), hence 

A = -^^ {wl'^" cos(0) + w/"^ sin(0) + + uoT^") = (19) 
2 uq 2 

Therefore, the Euler-Lagrange equations for the variables ^0', uq, u\\,'fl, r'^^ achieve the sim- 
ple expressions 

djl = ^ fi = const. (20) 
Tf,dr'^ - uods = ^ uo = -^V' (21) 
If.dr'f' + uiids = ^ U|| = — —l^, (22) 

d(j)' + qBds = ^ ^ = -fi, (23) 
ds 

^F'^Jr"" - l^du\\ - r^du^ - 2xqtiH,^ = 0, (24) 

ul-ul- 2qBti -1 = 0, (25) 

where VL = qB is the particle gyrofrequency (or cyclotron frequency) defined with respect 
to the rest mass. It is interesting to point out that Eq. (l23|) is formally identical to the 
non-relativistic one. Hence, at this order no correction due to space-time curvature appears 
in this equation. The remaining equations have a straightforward meaning. Thus, Eq. (|20l) 
implies the conservation of the relativistic magnetic moment /x, which is therefore an adiabatic 
invariant. Actually the result is expected to hold at any order in e and A, although the 
perturbative theory locally may not converge due to phase-space resonances. Furthermore, 
Eqs.f l23p and flMl) deliver the equations for the guiding-center 4- velocity, while fl25l) implies 
that the physical realizability condition is satisfied also by the guiding-center velocity. We 



stress that these equations, while equivalent to those obtained earlier |l|, have a simpler 
structure. This follows from the adoption of the synchronous variational principle here 
adopted which makes explicit the inner structure and relationship existing between the 
gyrokinetic variables. 



V. 5. CONCLUSIONS 



In this paper we have formulated a covariant gyrokinetic approach for single-particle dy- 
namics, which applies to curved space-time. The theory goes beyond the domain of validity 
of the treatment previously considered by Beklemishev and Tessarotto [l| and includes the 
case in which the electric field can locally vanish or result much smaller than the mag- 
netic field. As an illustration of the theory the perturbative theory has been developed to 
leading order in e. For this purpose we have adopted a synchronous variational principle 
based on constrained dynamics. The resulting Euler-Lagrange equations satisfy the required 
constraint demanded by the physical realizability condition for the 4-velocity and displays 
the inner relationships between the gyrokinetic variables. We have found that to leading- 
order no effect of curvature of space-time appears on the gyrofrequency at order e^. This 
information is relevant for plasma diagnostics in astrophysical plasmas. 
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